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I. INTRODUCTION 



The complex action of QCD at finite chemical potential makes it difficult to study finite- 
density QCD properties directly from first principle lattice calculations. There are many 
phenomenological models based on hadron degrees of freedom, such as the Walecka model 
[1], the Zimanyi-Mozkovski model [2] and the nonlinear a — uo model [3], as well as models 
which explicitly include quark degrees of freedom, for example, the quark meson coupling 
model [4], the cloudy bag model [5], the quark mean field model [6] and the NJL model [7]. 
With the development of these phenomenological models, the interactions between hadrons 
become more and more complex. The symmetries of QCD can be used to determine largely 
how the hadrons should interact with each other. With this in mind, the chiral SU(2) x 
SU(2) effective quark model was proposed. It has been used widely to investigate nuclear 
matter and finite nuclei at both zero and finite temperature [8]- [15]. 

In the last twenty years, exploring systems with strangeness, especially with a large 
strangeness fraction, has attracted a lot of interest. Therefore, a model which includes 
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strange quarks or hyperons is needed. Papazoglou et al. [16,17] proposed a chiral 577(3) 
model which extended the SU(2) chiral symmetry to S77(3)x 577(3). Recently, a chiral 
577(3) quark mean field model based on quark degrees of freedom was proposed by Wang 
et al. [18,19]. In this model, quarks are confined in the baryons by an effective poten- 
tial. The quark-meson interaction and meson self-interaction are based on SU (3) chiral 
symmetry. Through the mechanism of spontaneous symmetry breaking the resulting con- 
stituent quarks and mesons (except for the pseudoscalars) obtain masses. The introduction 
of an explicit symmetry breaking term in the meson self-interaction generates the masses 
of the pseudoscalar mesons which satisfy the relevant PCAC relations. The explicit sym- 
metry breaking term in the quark-meson interaction gives reasonable hyperon potentials in 
hadronic matter. This chiral SU (3) quark mean field model has been applied to investigate 
nuclear matter [20], strange hadronic matter [18], finite nuclei, hypernuclei [19], and quark 
matter [21]. By and large the results are in reasonable agreement with existing experimental 
data. 

A surprising result, which is quite different from most other models is that at some 
critical density the effective baryon mass drops to zero [22]. There is a first order phase 
transition around this critical density where the physical quantities change discontinuously. 
We will see later that this behavior is primarily a consequence of the nonlinear ansatz for 
the effective baryon mass which is related to the subtraction of the centre of mass (cm.) 
motion. In Ref. [23], the authors provided an exact solution of the effect of cm. motion and 
found that it was only very weakly dependent on the external field strength for the densities 
of interest. As a result, the cm. correction can be included in the zero point energy. In 
this paper, we will use this alternative definition of the effective baryon mass in medium to 
reexamine the properties of hadronic matter. 

The paper is organized as follows. The model is introduced in section II. In section III, we 
study infinite hadronic matter, finite nuclei, and hypernuclei with the improved treatment 
of the centre of mass correction to the baryon mass in medium. The numerical results are 
presented in section IV and section V is the summary. 



II. THE MODEL 

Our considerations are based on the chiral 577(3) quark mean field model (for details see 
Refs. [18,19]), which contains quarks and mesons as basic degrees of freedom. In the chiral 
limit, the quark field q can be split into left and right-handed parts q^ and q R . q = q^ + qn- 
Under 5£/(3)lX SU(3)r they transform as 

QL^q' L = Lq L , q R ^q' R = Rq R . (1) 

The spin-0 mesons are written in the compact form 

1 8 

M(M+) = E ± ill = -= (s a ± ip a ) A a , (2) 

v2 a=0 

where s a and p a are the nonets of scalar and pseudoscalar mesons, respectively, X a (a = 
1, 8) are the Gell-Mann matrices, and A = I. The alternatives, plus and minus signs 
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correspond to M and M + . Under chiral SU(3) transformations, M and M + transform as 
M - • \l' LMR+ and M+ -> M+ = RM+L+. The spin-1 mesons are arranged in a 
similar way as 

iM = \ (v, ± A.) = ^ E (*£ ± r ( 3 ) 

with the transformation properties: l^—^l'^ — Ll^L + , — > = Rr fM R + . The matrices E, 
II, and 71^ can be written in a form where the physical states are explicit. For the scalar 
and vector nonets, we have the expressions 



v2 a=0 



K* 



K* + \ 



*0 



if* 



(4) 



V * a =0 



(5) 



Pseudoscalar and pseudovector nonet mesons can be written in a similar fashion. 
The total effective Lagrangian is written: 



-eff 



£q0 + £aM + £sE + £-VV + £y-SB + £ 



J qM 



-XSB 



'Am £ 



(6) 



where £ g0 — Q il^d^ q is the free part for massless quarks. The quark-meson interaction C qM 
can be written in a chiral 577(3) invariant way as 



J-qM — 9s 



v/2 



/ 8 8 \ /8 8 \ 

* E + nTMJ * - 7 M E <A a - 7V E a?A a ) tt. (7) 

\a=0 a=0 / ^V^ \ a =0 a=0 / 



In the mean field approximation, the chiral-invariant scalar meson £ss and vector meson 
Cyv self-interaction terms are written as [18,19] 



£ee = ~\ koX 2 + C 2 ) + h (<x 2 + C 2 ) 2 + h + C 4 ) + h X a\ 



-hx 4 - 7X 4 In ^4 + -X 4 In 

4 Xo 

ix 2 



4 3 X ln a 2 Co' 



£ vv = 7 ^2 ( m '^ 2 + m > 2 + m l ( f )2 ) + ^ f^ 4 + 6cu 2 p 2 + p 4 + 20 4 ) , 
2 Xo v 7 v 7 



(8) 
(9) 



where <5 = 6/33; <7o, Co and xo are the vacuum expectation values of the corresponding mean 
fields a, ( and x- 

From the quark-meson interaction, the coupling constants between scalar mesons, vector 
mesons and quarks have the following relations: 
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9s 
9v 

2V2 



-q d 



- 9pP - 



'9p0 — 9u — 9u - 



1 s 
J_ Q s 
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ao 



g>=g» = g' = 0, 

9 s . = 9> =g; = gt = o. 



(10) 

(ii) 



Note, the values of <jq, (q and xo are determined from a minimization of the thermodynamic 
potential. On the other hand, the parameters a and Co are constrained by the spontaneous 
breaking of chiral symmetry and are expressed by the pion (F n = 93 MeV) and the kaon 
(F K = 115 MeV) leptonic decay constants as: 



Co = ^=(F n - 2F K ) 



(12) 



The Lagrangian C x sb generates the nonvanishing masses of pseudoscalar mesons 



p X 

L XSB ~ -2 
AO 



m\F v a + \V2m? K F K - ^LfA ( 



(13) 



leading to a nonvanishing divergence of the axial currents which in turn satisfy the partial 
conserved axial-vector current (PCAC) relations for re and K mesons. Pseudoscalar, scalar 
mesons and also the dilaton field x obtain mass terms by spontaneous breaking of chiral 
symmetry in the Lagrangian (8). The masses of u, d and s quarks are generated by the vac- 
uum expectation values of the two scalar mesons a and (. To obtain the correct constituent 
mass of the strange quark, an additional mass term has to be added: 



£a™ s = -Am s qSq 



(14) 



where S = | (i — XgV^j = diag(0, 0, 1) is the strangeness quark matrix. Based on these 
mechanisms, the quark constituent masses are finally given by 



m u = m d = --^=<7 and m s = -g s ( + Am s . 



(15) 



The parameters g s = 4.76 and Am s = 29 MeV are chosen to yield the constituent quark 
masses m q = 313 MeV and m s = 490 MeV. In order to obtain reasonable hyperon potentials 
in hadronic matter, we include an additional coupling between strange quarks and the scalar 
mesons a and ( [18]. This term is expressed as 



Ch = {ha + h 2 C) ss . 



(16) 



In the quark mean field model, quarks are confined in baryons by the Lagrangian C c = 
— ^ Xc^? (with Xc given in Eq. (17), below). The Dirac equation for a quark field under 
the additional influence of the meson mean fields is given by 



-ia ■ V + Xc(f) + firn* 



(17) 



where a = 7°7 , (3 = 7 , the subscripts % and j denote the quark i (i = u, d, s) in a baryon 
of type j (j = N, A, E, S) ; xdj') is a confinement potential, i.e. a static potential providing 
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confinement of quarks by meson mean-field configurations. The quark mass m* and energy 
e* are defined as 



m* = - g\C + m i0 (18) 

and 

e* = e t - glu - gty , (19) 

where is the energy of the quark under the influence of the meson mean fields. Here 
m i0 — for i — u, d (nonstrange quark) and m i0 = Am s = 29 MeV for i = s (strange 
quark). Using the solution of the Dirac equation (17) for the quark energy e* it has been 
common to define the effective mass of the baryon j through the ansatz: 



M* = y/Ef- < P f cm > , (20) 

where E* = J^i^e* + Ej spin is the baryon energy and < p*j 2 cm > is the subtraction of the 
contribution to the total energy associated with spurious center of mass motion. In the 
expression for the baryon energy riij is the number of quarks with flavor "i" in a baryon 
with flavor j, with j — N {p, n} , E {X^, S } , S E~} , A and Ej spin is the correction to 
the baryon energy which is determined from a fit to the data for baryon masses. 

There is an alternative way to remove the spurious cm. motion and determine the 
effective baryon masses. In Ref. [23], the removal of the spurious cm. motion for three 
quarks moving in a confining, relativistic oscillator potential was studied in some detail. It 
was found that when an external scalar potential was applied, the effective mass obtained 
from the interaction Lagrangian could be written as 

M; = ^n ije *-El (21) 

i 

where E® was found to be only very weakly dependent on the external field strength. We 
therefore use Eq. (21), with E® a constant, independent of the density, which is adjusted to 
give a best fit to the free baryon masses. 

Using the square root ansatz for the effective baryon mass, Eq. (20), the confining po- 
tential Xc is chosen as a combination of scalar (S) and scalar-vector (SV) potentials as in 
Ref. [19]: 

Xc(r) = i[XcW + X c SV W] (22) 



with 

-v s (r) = 

and 



X S c (r) = ]k c r\ (23) 



Xc SV (r) = ^ c r 2 (l + 7 °). (24) 

On the other hand, using the linear definition of effective baryon mass, Eq. (21), the confining 
potential Xc is chosen to be the purely scalar potential Xc( r )- The coupling k c is taken as 
k c — 1 (GeV fur 2 ), which yields baryon radii of about 0.6 fm. In both cases, the baryon 
masses in vacuum are chosen as: Mn = 939 MeV, Ma = 1116 MeV, Ms = 1196 MeV and 
M= = 1318 MeV. 
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III. HADRONIC SYSTEM 



F 


= 9^uj u 










E^ 


= d,A u 


- d^. 



Based on the previously defined quark mean field model the effective Lagrangian for the 
study of hadronic systems is written as 

C H = Mi^d, - M B )^ B + ^d»od»o + ^ M C^C + \d, X d»X - - 

—E^E^ - g^Bi^B^ - g^Bi^B^ - \g B p i>Bi^Bff 

-e B ^Bl^ B A^ + C M , (25) 

where 

(26) 
(27) 
(28) 

tB is the coupling constant of comloub interaction. The mesonic Lagrangian 

Cm = £ee + Cyy + C x sb (29) 

describes the interaction between mesons which includes the scalar meson self-interaction 
£ee, the vector meson self- interaction Cyy an d the explicit chiral symmetry breaking term 
C x sb defined previously in Eqs. (8), (9) and (13). The Lagrangian Cm involves scalar 
(a, ( and x) an d vector (u and 0) mesons. The interactions between quarks and scalar 
mesons result in the effective baryon masses M B , where subscript B labels the baryon flavor 
B = N, A, E or S. The interactions between quarks and vector mesons generate the baryon- 
vector meson interaction terms. The corresponding vector coupling constants g% and gjjf 
satisfy the 377(3) flavor symmetry relations: 

A E o S ^ N u g v A E ^ 3 N s 9v /orA 

9 U =9 U = 2 9u = 3^ = 2 9u, = "y= and 9^=9^ = ^9<p = -^9 U =9^ = ^- (30) 

Compared with the chiral 577(3) model of Refs. [16] and [17], the mesonic Lagrangian is the 
same. However, that model is based on the hadron degree of freedom and the effective baryon 
masses are obtained from the direct baryon-meson coupling. This chiral £77(3) quark mean 
field model is on the quark level and the effective baryon masses are generated by quarks 
which couple with mesons and are confined in baryons by a potential. In the earlier work, 
the effective baryon masses are calculated with the square root ansatz. In this paper, we 
use the improved linear definition of baryon masses. 

The equations for mesons 0j can be obtained by the formula — d ^ 9 H ^ — ^f- = 0. There- 
fore, the equations for a, ( and x are 

-d^a + fcoxV - 4fci (a 2 + ( 2 W - 2k 2 a 3 - 2k 3X a( - + —^F^ 

3a Xo 



— m w uj — h > — — < ip B iPB >= 0, (31) 
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-d^C, + k oX 2 ( - 4fcx (a 2 + C 2 ) C - 4A; 2 C 3 - W - A/ + ^ (^m 2 ^ - ^m*^) 



- ( — I m^0 2 
VW 



(32) 



-d^ X + koX (a 3 + C 2 ) " h° 2 ( + ( 4A; 4 + 1 + 4 In 



+ 



2x 

xl 



x_ 

Xo 

mlF n a + (V2m 2 k F k - 7^™^) C 



--ln^A 

X 2 2 
2 m . W : 

Xo 



= 0. 



(33) 



For infinite hadronic matter, the meson mean field is independent of position and < ips^B > 
is expressed as 



9b 



[ Mb r 
it 2 Jo 



dk 



M* B 2 + k 2 



(34) 



9b M , 
2tt 2 



*3 r 



M* B \ 



1 + 



k 2 
M* B 2 



In 



K F B 

Ml 



+ 



1 



k F B 

M* 2 



For finite nuclei, the integration in < ipsipB > will change to a sum over the states ip%, 
which are obtained from the Dirac equation for the baryon B in state a 

1 



-«7 • V + M* B + u<y + g^ 07 O + -g p P7o + e B ^4o7o 



(35) 



This set of coupled equations is solved iteratively. At each iteration, we first solve the 
equations for the quark wave function using fourth-order Runge-Kutta for given meson fields 
and obtain the effective baryon masses. We then solve the equations for the nucleon radial 
wave functions with the same method. The corresponding scalar and vector densities can 
be obtained. Then the mean meson fields which will be used in the next iteration step 
are determined for the calculated densities. The iteration is stopped when convergence is 
achieved. The energy for the finite system within the mean field approximation can be 
derived in the standard way. 

n 1 - dM* dM* 



E 



Si ! r dM* dM* 

£ _ 4(2*. + 1) - 2 j dM(°-^- + c^)pi + M + q 3 M + 9 j M) + 



j=N,A,~;a=l 



E r 



(36) 



where fP s = ipjipj and pj = ^70^ ■ ^ a are the Dirac single particle energies and j a are the 
total angular momenta of the single particle states. By using the equations of meson fields, 
the rearrangement energy E rearr can be written as 



E rearr = J d 3 r j^cu 4 + p 4 + 6cu 2 p 2 + 20 4 ) - h(a 2 + ( 2 ) 2 - k 2 ^ 



a 



■l^hxv C + 3X In 



^2^2 / + 9 ^2 m TtJK& 

°oso / z Xo 



(^ m ^^ _ % 7f m ^^) VvaC) 



(37) 
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where the constant V vac is the vacuum energy which is subtracted to yield zero energy in 
the vacuum. 



IV. NUMERICAL RESULTS 

Before doing detailed numerical studies, we first determine the parameters in the model. 
The coupling constant g s as well as Am s0 are determined by the constituent quark masses, 
m q and m s . m v and p, are obtained by fitting the vector meson masses. The confining 
coefficient, k C) is chosen to be 1000 MeVfm -2 to make the baryon radii (in the absence of a 
pion cloud [24]) around 0.6 fm. For nuclear matter, there are seven other parameters, ko, k±, 
k 2 , &3, £4, #4 and g v , to be determined. We fit them to the 7r-meson mass, A"-meson mass 
and the average mass of rj and rj' which are given by the eigenvalues of the mass matrix 



5 2 C 



H 



6<f>i6<f>j 



(38) 



There are two constraints associated with the saturation properties of nuclear matter. We 
choose the parameters to fit the binding energy of nuclear matter e/p — = —16 MeV 
at saturation density po — 0.16 fm -3 . The parameters should also produce a reasonable 
compression modulus and effective nucleon mass at saturation density. For strange matter, 
there are two additional parameters hi and hi to be determined. They are restricted by 
the hyperon potentials in hadronic matter. All these nine parameters are listed in Table I 
with the two methods for computing the effective baryon mass. The properties of nuclear 
matter and the masses of a and ( are listed in the Table II. From the table, one can see that 
the main difference between the two versions is that in the linear definition (Eq. (21)), the 
effective mass is larger. In fact, if we choose a larger value of g±, the effective mass will be 
smaller. However, in that case the results for finite nuclei will not be so good. 

TABLE I. Parameters of the model. 



version 


k 


ki 


k 2 


h 


/C4 


9s 


9v 


9a 


hi 


h 2 


square 


4.21 


2.26 


-10.16 


-4.38 


-0.13 


4.76 


10.99 


7.5 


-2.07 


2.90 


linear 


3.97 


2.18 


-10.16 


-4.15 


-0.14 


4.76 


8.70 


15.0 


-2.66 


2.45 



TABLE II. Nuclear properties and scalar meson masses. 



version 


po (fm" 3 ) 


E/A (MeV) 


M* N /M N 


K (MeV) 


m a (MeV) 


m c (MeV) 


square 


0.16 


-16.0 


0.603 


225 


466.2 


1167.1 


linear 


0.16 


-16.0 


0.742 


303 


487.8 


1168.0 
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A. Infinite hadronic matter 



We first consider infinite nuclear matter. In Fig. 1, we plot the effective nucleon mass 
versus meson mean field a. The solid and dashed lines correspond to the linear (21) and 
square root (20) definitions of baryon mass, respectively. In vacuum, we have <7 = F n ~ 0.47 
fm _1 . With increasing density the value of a increases, resulting in a decreasing effective 
nucleon mass. For a > —0.2 fm _1 , the dashed line decreases very fast, while the slope of the 
solid line changes a little. In Fig. 2, we show the effective nucleon mass divided by the free 
nucleon mass, M^/M N , versus nuclear density. The effective mass for the square root ansatz 
decreases faster than that for the linear definition and, as a consequence, at some critical 
density the effective mass drops to zero. However, in the case of the linear definition, the 
effective mass decreases slowly at high density and there is no phase transition to a state of 
chiral symmetry restoration. We plot the energy per nucleon versus nuclear density in Fig. 
3. The behavior of the solid line and dashed lines are close when the density is small, say 
p B < 0.2 fm~ 3 . Both curves pass through the saturation point of nuclear matter. For the 
square root ansatz for effective mass, E/A changes disco ntinuously at the critical density. 

n i 1 i 1 i 1 i 1 i 1 i 1 i I 




FIG. 1. The effective nucleon mass M^/Mn versus the a mean field The solid and dashed lines 
are for linear and square root treatments of the effective baryon mass, respectively. 
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0.1 0.2 0.3 0.4 0.5 0.6 
Pn ( frrf 3 ) 

FIG. 2. The effective nucleoli mass M^/Mjv versus nuclear density pj^. The solid and dashed 
lines are for linear and square root treatments of effective baryon mass, respectively. 
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FIG. 3. The energy per nucleon E/A versus nuclear density pN. The solid and dashed lines are 
for linear and square root treatments of effective baryon mass, respectively. 



TABLE III. Hyperon potentials in MeV. 



version 


U {N) 

N 


rr(N) 
U A 




ui N) 


rr(A) 




c4 A > 


MS) 

U A 




u<P 


square 
linear 


-74.1 
-64.0 


-28.5 
-28.0 


-22.8 
-28.0 


-14.3 
8.0 


-31.7 
-24.5 


-29.3 
-24.5 


-33.2 
-20.6 


-41.8 
-30.6 


-36.9 
-30.6 


-40.1 
-50.3 
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0.1 0.2 0.3 0.4 0.5 0.6 
Pb ( fm" 3 ) 

FIG. 4. The effective baryon mass versus baryon density ps calculated for the linear definition 
of baryon mass with strangeness fraction f s = 0. The solid, dashed, dotted and dash-dotted lines 
are for nucleon, A, £ and S, respectively. 

Now we consider strange hadronic matter which includes A, E and S hyperons. Since 
we have studied strange hadronic matter with the square root ansatz for the baryon mass 
in previous work [18], we now concentrate on the linear definition and compare the results 
of these two cases. Before studying strange matter, one should first reproduce the hyperon 
potentials in hadronic matter. The hyperon potential felt by baryon j in i-matter is defined 
as 

= M* - M, + gico + g^. (39) 

In Table III, we list the hyperon potentials. For the A hyperon, the empirical value of 
at the saturation density of nuclear matter, p , is —28 MeV [25,26]. For recent 
experiments suggest that may be —14 or less [27,28]. In A matter, the typical values of 
L7( A ) (j = A, S) are around -20 MeV at density p = p /2 [29]. In S matter, (j = A, S) 
are around —40 MeV at density p = p [29]. From Table III, one can see that in the 
case of the square root ansatz for baryon mass, the S potential in hadronic matter is more 
reasonable, while for the case of linear definition, the A potential is more reasonable. For 
the S potential in nuclear matter, the results are consistent with the earlier analysis [30]. 
Calculations in Brueckner-Hartree-Fock method show that the potential in symmetric 
nuclear matter is about -20 MeV [31]. However, the modern fits to £ atomic data suggest a 
repulsive potential [32]. Recent experiment also shows that a strongly repulsive S-nucleus 
potential is required to reproduce the shape of the (n~ , K + ) spectra on nuclear targets 
[33]. From the equation (39), one can see that the hyperon potentials are determined by the 
decrease of the effective baryon mass and the vector interaction. Since A and £ have the 
same quark components, the calculated potentials and in this model are close. 
It is important and interesting to understand why there is a large difference of A and £ 
potentials which eventually could be considered in the model. The precise value of will 
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not change the main results of the paper because S hyperons do not appear in the strange 
hadronic matter untill very large baryon density. This is because the chemical potentials of 
A and E are the same, while the mass of E is 80 MeV larger than that of A. 
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Pe ( fm' 3 ) 

FIG. 5. The effective baryon mass versus baryon density ps calculated for the linear definition 
of baryon mass with strangeness fraction f s = 2. The solid, dashed, dotted and dash-dotted lines 
are for nucleon, A, £ and S, respectively. 
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FIG. 6. The energy per nucleon E/A versus baryon density pb calculated for the linear definition 
of baryon mass with differen strangeness fraction f s . 



In Fig. 4, we plot the effective baryon masses versus density with strangeness fraction 
f s = with the linear definition of mass. f s is defined as 
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PB 

where pa, ps, p~ are the baryon densities of A, E and S, respectively, and p# is the total 
density of all kinds of baryons. All the baryon masses decrease smoothly with increasing 
baryon density. The nucleon mass drops faster than other baryons and the mass of S drops 
slowly. This is because in nuclear matter, the interaction between nonstrange quarks is 
stronger than with the strange quark. With increasing strangeness fraction, the interaction 
between strange quarks becomes more important. We show in Fig. 5 the effective baryon 
mass versus density with f s = 2. In contrast with Fig. 4, the mass of the S hyperon drops 
faster than that of other baryons since the S has more strange quarks. For any f s , the baryon 
masses decrease slowly and smoothly at high density which is different from the case of the 
square root ansatz for baryon mass, where the baryon mass changes discontinuously at some 
high density. The energy per baryon versus density with different strangeness fractions is 
shown in Fig. 6 for the linear definition of baryon mass. With increasing f s , the binding 
energy of strange hadronic matter increases. The maximum binding energy is about 108 
MeV, where the corresponding f s is about 1.97. For the square root ansatz for baryon mass, 
the largest binding energy is about 70 MeV and the corresponding f s is about 1.5 [34]. 
The results are comparable with those of model N in Ref. [35] where the maximum binding 
energy is about 79 MeV with f s = 1.45. The binding energies of various { N, A, S} and 
{N, A, E, S} systems were also calculated using the Brueckner-Hartree-Fock approximation 
[36]. Our results are close to the case of Ref. [36] where the ratios of A, E and S in hadronic 
matter are the same. 



B. Finite nuclei and hypernuclei 



We now investigate the finite system. We do not adjust the parameters, rather they are 
the same as in infinite hadronic matter. The charge density versus the radius is shown in 
Fig. 7. Both of the effective baryon mass definitions give reasonable results through there 
is a little difference from the experiment when the radius is smaller than 2 fm. In Fig. 8. 
the charge density of 208 Pb is shown. Again, the two definitions give similar results. This 
arises naturally since the nucleon density is around saturation density in the center of the 
finite nuclei and both treatments yield the correct saturation properties of nuclear matter. 
We plot the proton energy levels of 16 and 208 Pb in Fig. 9 and Fig. 10. The energy levels 
are qualitatively reproduced. For the square root ansatz for baryon mass, the spin-orbit 
splitting is quite close to the experiments. For example, the proton spin-orbit splitting of 
lpi/2 and lps/2 of 16 is about 5.5 MeV which is close to the experimental value 6.7 MeV. 
The splitting of lg 7/2 and lg 9/2 ( 2d 3/2 and 2d 5/2 ) of 208 Pb is about 4.2 MeV ( 1.6 MeV ) 
which is close to the experimental value 3.9 MeV ( 1.5 MeV ). For the linear definition of 
baryon mass, the spin-orbit splitting is smaller. This is because at saturation density, the 
effective nucleon mass is higher in this definition and the spin-orbit splitting is proportional 
to the decrease of the effective nucleon mass. However, the smaller spin-orbit splitting can 
be improved by going beyond the mean field approximation [37]. 
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FIG. 7. The charge density versus radius for 16 O. The dashed and dotted lines are calculated 
with linear and square root treatments of effective baryon mass, respectively. The solid line is from 
the experimental data. 
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the experimental data. 
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FIG. 9. The proton energy levels for 16 O. The first and second columns are calculated with 
linear and square root treatments of effective baryon mass, respectively. The third column is from 
the experimental data. 
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FIG. 10. The proton energy levels for 208 Pb. The first and second columns are calculated with 
linear and square root definition of effective baryon mass, respectively. The third column is from 
the experimental data. 
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Next we investigate hypernuclei. The hypernuclei which includes one hyperon was stud- 
ied in some detail in the QMC model [38] . We will concentrate on the binding energies of A 
and double-A hypernuclei. The binding energy B\, of A-hypernuclei is expressed as 

B A = M( A - 1 Z) + M A -M(iZ). (41) 

The results are plotted in Fig. 11. The experimental values are from Refs. [39,40]. For the 
square root ansatz of baryon mass, one can see that for the light hypernuclei, the binding 
energies B A are about 3 MeV larger than the experimental values. When the baryon number 
is larger than 10, the deviation from the experimental values is around 20-30%. For the 
heavy lambda-hypernuclei, the results are very close to the experiment. This is because 
the parameters are obtained for bulk hadronic matter and the mean field approximation is 
not good when baryon number A is not large. For the linear definition of baryon mass, the 
results are improved since the A potentials are more reasonable in this definition. 
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FIG. 11. Hypernuclei binding energies B\ calculated with linear and square root treatments of 
effective baryon mass versus baryon number A. 



TABLE IV. The values of double-A hypernuclei B AA and AB\ A in MeV. 



nuclei 


5aa(L) 


Saa(S) 


Exp. [41] 


Exp. [42] 


A£? AA (L) 


A£?aa(S) 


Exp. [41] 


Exp. [42] 


AAHe 


11.45 


15.21 


10.9±0.8 


7.25±0.19 


1.61 


1.97 


4.7±1.0 


1.01±0.20 


A° A Be 


19.66 


23.66 


17.7±0.4 




1.32 


1.60 


4.3±0.4 




aaB 


23.67 


30.22 


27.5±0.7 




0.89 


0.54 


4.8±0.7 
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The binding energies of two lambdas _B AA defined as 

B AA = M( A - 2 Z) + 2M A - M(i A Z) (42) 

are listed in table IV. One can see that for the square root ansatz for baryon mass, the 
calculated results are several MeV larger compared with the old experimental values [41]. 
The linear definition improves this result. There are two reports for the new events of the 
production and detection of the double lambda hypernuclei [42,43]. It shows that B AA of 
A 6 A He is much smaller than the old value [42]. This new event was discussed in detail in Ref. 
[44]. The A-A interaction energy AB AA is defined as AB AA = B AA — 2B A . The calculated 
results of these two treatments are similar and much smaller than the old experimental 
values. AB AA of A 6 A He obtained in this model is comparable with the new result. 



V. SUMMARY 

We have used an improved treatment of the cm. motion in calculating the effective, 
in-medium, baryon mass in an investigation of infinite hadronic matter, finite nuclei and 
hypernuclei within the chiral SU(3) quark mean field model. The results are compared 
with earlier results which used the square root ansatz for effective mass. Both treatments 
fit the saturation properties of nuclear matter and therefore, for densities lower than the 
saturation density, these two treatments give reasonably similar results. The S potential in 
hadronic matter is reproduced better in the square root ansatz for baryon mass, while the 
linear definition gives a better A potential. As a result, the binding energy of A hypernuclei 
calculated from the linear definition is better when compared with experimental values. 
The energy levels of finite nuclei are reasonable in both of the treatments. In the linear 
definition, the spin-orbit splitting is smaller than that in the square root case. This is 
caused by the different effective nucleon mass at saturation density. The spin-orbit splitting 
can be improved by going beyond the mean field approximation [37]. 

For high baryon density, the predictions of these two treatments are quite different. There 
is a phase transition of chiral symmetry restoration in the case of the square root ansatz for 
baryon mass. The physical quantities, such as effective baryon mass and energy per baryon 
change discontinuously at the critical density. In the linear case, no such transition occurs. 
The effective baryon masses decrease slowly at high density. The different behavior at high 
density will result in significant difference for high density physics. It is therefore of interest 
to study the properties of neutron stars with these two treatments and to try to construct 
clearer theoretical motivations for possible definitions of the effective baryon mass. This will 
be done in the future. 
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